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Although both of the subject theories have shortcomings, they have been shown to give 
meaningful results in a variety of applications. For instance, low-aspect-ratio theory has been 
shown to give reasonable results for the steady state aerodynamic characteristics of slender config- 
urations at transonic speeds, 6 and for determining the static divergence characteristics of all move- 
able, delta-wing control surfaces. 7 Piston theory has been used in flutter calculations of a variety 
of configurations. 8 * 9 Indeed, piston theory has been routinely used in parametric aeroelastic analy- 
ses of advanced configurations, most recently in National Aero-Space Plane (NASP) applica- 
tions. 10 Even in instances where quantitative results may be inaccurate, trend data obtained by 
using these theories are usually correct. 

The present author believes that there is much to be learned by just going through the pro- 
cess of using aerodynamic sensitivity derivatives in multidisciplinary optimization. The knowledge 
gained from formulating and integrating aerodynamic sensitivity derivatives into the multidisci- 
plinary process has significant value in and of itself. It is not necessary that the aerodynamic the- 
ory used be of "CFD quality", it is only necessary that the aerodynamic theory give reasonable re- 
sults for the problem being addressed. Of course, in the long term, accurate CFD methods are 
needed before the full benefits of multidisciplinary optimizations can be realized in the design of 
future aerospace vehicles. Sobieski's plea is, therefore, still relevant and hopefully being heeded 
by CFD researchers. Nevertheless, it is hoped that the present exposition will make it clear that it 
is not necessary to await further aerodynamic theory developments before proceeding with the use 
of aerodynamic sensitivity derivatives in multidisciplinary optimization research. 

AERODYNAMIC THEORIES 

Two existing aerodynamic theories for which it is relatively simple and straightforward to 


Low-Aspect-Ratio Theory 

Planar wing- steady flow : - Low-aspect-ratio aerodynamic theory as used herein is appli- 
cable to planar, slender, pointed, symmetric wings such as that shown in Sketch 1. The shape of 
the wing is defined by the function S(x) which is monotonically increasing in the stream direction. 
This requirement is necessary because portions of a wing downstream of the point of maximum 
span do not generate any lift. The root chord is denoted by C r , and the semispan at the trailing 
edge by Sq. The wing is inclined at an angle a to the free stream velocity V. The flow is consid- 
ered to be two dimensional transverse to the flow direction. Compressibility effects are considered 
to be small. 
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Sketch 1 


As developed by Jones 3 , the lifting pressure distribution in the spanwise direction y at 
streamwise station x for a wing at angle of attack a is given by 

p(x,y) - 2 a pV2 = ^T « « P v2 ^[VS(x)2-y 2 ]. (1) 

The lift per unit length at streamwise station x is obtained by integrating equation 1 across the span. 
The result is 

^ = 2k a p V2 S(x) ^ = jc a pV2 ^ [s2(x)] . (2) 
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The familiar expression for the total lift on the slender pointed wing results from the 
streamwise integration of equation 2 from the wing apex, x=0, to the trailing edge, x=Cf. The re- 
sulting expression for the total lift is 


C r 

L = n a pV 2 J ^ [s 2 (x)] dx = n a pV 2 S 0 2 . (3) 

0 

The total lift depends only on the width (maximum span) of the wing and not on the area, nor does 
it depend on the detailed shape of the planform. 

The aerodynamic moment about the leading edge may be obtained by multiplying equa- 
tion 2 by x and then integrating the result over the chord from the wing apex to the trailing edge. 
The resulting expression for the moment is 


C r 

M = 7t a pV 2 { CrS 0 2 - J S 2 (x) dx } . (4) 


Equations 3 and 4 are the basic equations defining the steady aerodynamic lift and moment on a 
low-aspect-ratio wing of general planform. 

These lift and moment equations can be made more planform specific, yet still retain quite a 
bit of generality, if the planform shape function is assumed to be of the form 

S(x) = K x n (5) 

where K is a constant and n is an integer number. Substituting this shape function into equations 3 
and 4, respectively, yields for the lift 

L = k a pV 2 K 2 C r 2n , (6) 


and for the moment 


M = n a pV 2 K 2 ) C r 2n+l - L (jn+j ) Cr- 


(7) 



For the special case of a triangular wing, K = tanO (The angle 0 is the complementary angle 
to the sweep angle.) and n=l, the expressions for the lift and moment become, respectfully, 


L = 7t a pV 2 tan 2 © Q 2 , 


( 8 ) 


and 


M = ^ n a pV 2 tan 2 0 Q 3 = L ( j C r ) (9) 

It is readily seen that the center of pressure is at two-thirds of the root chord measured from the 
apex. 

An interesting result is obtained if the lift for the triangular wing is converted to coefficient 
form by dividing equation 8 by the product of the dynamic pressure ^ pV 2 and the wing area 

C r S 0 . The resulting expression for the lift coefficient is 

Cl = |aRoc (10) 

4S 2 

where the wing aspect ratio AR =of^- The pitching moment coefficient can be obtained by divid- 
ing equation 9 for the moment by the dynamic pressure, wing area, and reference chord C r . The 
resulting expression for the pitching moment coefficient is 

C M =f ARa = (|) C L . (ID 

It is possible to develop aerodynamic sensitivity derivatives that are functions of selected 
planform variables by differentiating the appropriate equation(s) presented above. For example, 
the variation of total lift with maximum semispan (semispan at the trailing edge) can be obtained by 
differentiating equation 3. The result is 

^ = 27t a pV 2 S Q . (12) 
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Another example is to use equation 10 to determine the sensitivity of the lift coefficient to changes 
in aspect ratio. 

Deformed wing, unsteady flow: - Garrick 1 1 extended low-aspect-ratio theory to the un- 
steady case of a deforming (nonplanar) wing. In discussing this extension, consider Sketch 2 that 
illustrates a slender wing that is deformed by bending out of the plane of the wing. The deformed 
shape is denoted by Z(x,y,t). At a given chordwise station, the deformation is constant over the 
span. The normal wash w(x,y,t) at a point x, y on the surface of a wing at an instance of time t is 
given by the relationship. 
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The spanwise pressure distribution at a particular chordwise station as derived by Garrick 
is given by 

P<*,,> = *2p [ ^2 + + ] VW 


dx^ 


r 9Z(x,y,t) v 3Z(x,y,t) i S(x) dS(x) 

2pv[- a,— +v — 5, J q^rar- 


(14) 


Integration of this pressure distribution over the span gives the following result for the chordwise 
lift distribution. 
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%^=-xpS 2 (x)[^^ + 2V^^ + V 2 ^%^] 


3x 2 


9t3x 


3x 2 


- 2n pV S(x) [ aZ( ^’ y,t) + V~ Z ^ y,t) ]. 


3x 


( 15 ) 


Both the pressure and lift are functions of geometry, S(x) and the bending slope 

and curvature — the time rate of change of the bending slope and the 

^7/x y t) 0Zfx V 0 

acceleration — and velocity — y ’ of the bending deformations. 

For the steady state case (Z not a function of time t) the pressure distribution is given by 

p(x,y,t) = -2 pV 2 Vs(x) 2 -y 2 . 2 pV 2 3Z ^ X ' y) - S(x) S§&! 

3x2 VswV dx 


= -2pV^[^^W], 

and the chordwise lift distribution is given by 

= -7i pV 2 S 2 (x) a2Z ^’ y ’ t} - 2 tc pV 2 S(x) ^ 7 ^ dZ(x,y,t) 


(16) 


dx 3x 


= -,pv 2 3 [^S 2 W ]. 


(17) 


The total lift is obtained by integrating equation 17 from the leading edge apex to the trailing edge. 
The result is given by equation 18. 


L = -jt pV 2 S G 2 — ( ^ r,Q) . 


(18) 


Equation 18 is very similar to the expression for lift of the undeformed planar wing, equation 3, 
except now the angle of attack is the bending slope at the trailing edge. Thus, the total lift on a 
statically deformed low-aspect-ratio wing does not depend on the chordwise distribution of the de- 
formed shape, but rather, only depends on the slope at the trailing edge of the wing. 
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The expression for the moment about the leading edge is given by 


Q 


M = -7t pV 2 { CrSo 2 8Z( ^ r, ° ' - J 3Z( 3x y>t ’ S2(x) dx I* 

0 


( 19 ) 


This equation is very similar to the expression for the moment for the planar wing, equation 4. 

The lift and moment for a deformed wing that is at an initial angle of attack a are simply the 
respective sums of equations 3 and 18 and of equations 4 and 19. 

Aerodynamic sensitivity derivatives for the deformed low-aspect-ratio wing can be deter- 
mined from the equations presented in this section in the same manner as pointed out for the planar 
wing case. For example, the change in chordwise lift distribution with span can be obtained from 
differentiating equation 15. It should be noted that in practical optimization problems the 
displacement related terms in equation 15 may be functions of the span and this fact must be taken 
into account in performing the required differentiations to determine the aerodynamic sensitivity 
derivatives. By way of example, consider equation 18 for the total steady lift which is a simpler 
expression than equation 15 but still can be used to illustrate the point. If the displacement is not a 
function of the span, then the change in total lift with respect to the maximum span is given by 

BL . nV? 0 az(c r ,o) 

5 3- = -27tpV^S 0 — . (20) 


Now, if the displacement is a function of the span, the sensitivity becomes 

3L , , , 0 „ o 3Z(C r ,0) 3 2 Z(C r ,0) 

5s; - = - 2k p v2 V ,v - K p v s o asJK • 




(21) 


The first term in equation 21 is the same as equation 20. The second term results from the fact that 
the displacement is a function of the span. 


Piston Theory 

Piston theory 4 - 5 is a relatively simple supersonic aerodynamic theory. It has most often 
been used in supersonic flutter analyses, but it has steady state applications as well. Piston theory 
gets its name from the notion that the pressure at a point on an airfoil is related to the normal veloc- 
ity of the airfoil in the same way that the pressure on the face of a piston moving in a channel is 
related to the motion of the piston. There are limitations on the use of piston theory. It is consid- 
ered to be accurate when M 2 »l and when the product M8 is small. The term 8 is the larger of the 
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maximum thickness- to-chord ratio of the airfoil section or the ratio of the dynamic displacement to 
airfoil section chord. 

The relationship between the pressure p on the face of the piston and the motion of the 
piston w(t) in a channel as illustrated in Sketch 3 is given by 



Sketch 3 

where p Q is the pressure in the undisturbed fluid, a is the speed of sound in the fluid, and y is the 
ratio of specific heats. The expression given by equation 22 is the third-order relationship. Typi- 
cally only second-order piston theory is used in aerodynamic applications; the cubic term in equa- 
tion 22 is not used. Thus, the second-order expression for the pressure becomes 

PPo = pa 2 {[^] + ^Pf] 2 }- 03) 

In aerodynamic applications, the piston motion term w(t) in equation 23 is replaced by the 
normal wash of the airfoil. The normal wash on the upper surface of an airfoil surface is given by 

w(x,y,t) = + (Vj| + 1) Z(x,y.t) + V (24) 

where Z(x,y,t) is the position of the mean surface of the airfoil and z(x,y,t) is the function that de- 
scribes the airfoil contour measured from the mean surface. These terms are illustrated in Sketch 4. 
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The corresponding normal wash on the lower surface is given by 



Z(x,y,t) + V 


dz(x,y) 

dx 


( 25 ) 



Sketch 4 


The normal wash as given by equations 24 and 25 can be respectively substituted into 
equation 23 to provide expressions for the pressure on the upper surface and the lower surface of 
the airfoil. In most instances, however, it is the lifting pressure (difference between upper surface 
pressure and lower surface pressure) that is desired. By performing the two substitutions just 
mentioned and then taking the difference between the results, the lifting pressure is obtained. The 
resulting expression for the lifting pressure is 

Ap(x,y,t) = -2 pa [ 1 + T ^ z(x,y) ] [ (V^ + Z(x,y,t) ] (26) 


where T = M . The pressure at a point x,y is only dependent on the geometric characteristics 

of the airfoil section at that point and on the characteristics of the deformation at that point There 
is no mutual interaction of a disturbance at one point on the pressure at another point. 


It should be pointed out that there is an alternate form for equation 23 that follows from 
quasi-steady second-order linear supersonic theory. 12 This alternate form is 

p-po -.*<*) Wl* [ ^gr^ l p?] 2 } (27 > 

where P = Vm 2 - 1 . At high Mach numbers where P=M, this result approaches the second- 
order piston theory pressure given by equation 23. It is expected that the pressures given by equa- 
tion 27 will be more accurate at lower supersonic Mach numbers than the pressure given by equa- 
tion 23. The lifting pressure can be obtained by using the normal wash expressions given by 
equations 22 and 23 with the pressure equation 27. The resulting expression for the lifting pres- 
sure is 


Ap(x,y,t) = -2pa [ 1 + T ^ z(x,y) ][ 0^“ + Z(x,y,t) ] 


(28) 


but where now T = 


(Y+1)M 4 - 4p 2 
2p3 


. The form of this equation and equation 26 are the same. 


the only difference being in the value of the constant T. 

The proceeding equations are put down for the unsteady case, but can be readily converted 
to the steady flow case by deleting the time dependent terms. That is, the deformation shape Z in 
equation 26 is not a function of time. The expression for the pressure in steady flow is 


A p(x,y) = .2pa[l + r^i][v^] 


dx 


(29) 


where V is one of the two values listed previously, depending on whether the second-order piston 
theory form or the second-order supersonic theory form is desired. 

Because the surface pressures predicted by piston theory are a function of the airfoil sec- 
tion, this theory can be used to determine aerodynamic sensitivity derivatives that are functions of 
the airfoil section, for example, airfoil thickness-to-chord ratio. By way of illustration consider the 
simple delta airfoil illustrated in Sketch 5. The airfoil section is defined by the relationship 

z(x,y) = ( ^ ) 


X. 


(30) 



The slope of the airfoil is given by 



( 31 ) 



After the value of the airfoil section slope given by equation 31 is substituted into equation 29 and 
the derivative with respect to the airfoil thickness is taken, the sensitivity of the pressure at point 
x,y to the maximum thickness is determined. This expression is 


*Wp.-2pa[l + r(i)][(v 


3Z(x,y) \ 1 
dx > J 


(32) 


Sensitivity of forces and moments with respect to other geometric variables such as area 
can also be obtained in a straightforward manner. To illustrate this, consider the clipped delta wing 
at constant angle of attack a that is illustrated in Sketch 6. The steady state lift for this configura- 
tion can be determined for the steady flow case by integrating the pressure given by equation 26 
over the area of the planform, thus, 


Span Chord. 

L = -2 pa j J 

0 0 


[ |+r 3zfes) ][ v ^M)]d xd y. (33) 


dx 
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For this example the airfoil thickness-to-chord ratio is a constant independent of span, so 


^ = (i)=x. 


(34) 


Because the undeformed wing is at constant angle of attack a. 


^x,y) = a 


(35) 


By using the relationships in equations 34 and 35 with the lift expression, equation 33, then the 
total lift becomes 


S 

L = 2 paV a j 

0 


C r 

J 

X LE 


[ 1 + r t ] dx dy 


1 3 


(36) 



where x LE = ( Cr g C l ) y for the trapezoidal wing shown in the sketch. After the integration is 
performed and terms collected, the total lift is given by 

L-2 paV a [ 1 + T x ] Area 

= ^ pV 2 ^ a [ 1 + T x ] Area (37) 

where the ratio ^ containing the Mach number M is the familiar linear supersonic lift curve slope. 

The sensitivity of the total lift with planform area for the trapezoidal wing used in this example is 
simply the derivative of equation 37 with respect to the area, or 

55ta = I p V 2 Si “1 1 + rt l • O® 


The sensitivity of aerodynamic quantities with respect to other geometric parameters may be de- 
termined by similar means. As was pointed out for low-aspect-ratio theory case, the deformation 
shape, either static or dynamic, may be a function of the geometric parameter so this fact must be 
taken into account in performing the required differentiations. 

CANDIDATE PROBLEMS 

As mentioned previously, there are classes of problems for which low-aspect-ratio theory 
and piston theory would be expected to given reasonably accurate results. In this section, a couple 
of optimization problems for their use are suggested. Of the two theories, piston theory is the most 
attractive for use because it is applicable to a greater variety of planforms and more suitable to be 
used in combined steady and unsteady applications. 

The first problem is the design of a low-aspect-ratio all-moveable control fin that might be 
used on a missile. Both theories would be applicable to this problem. The objective would be to 
design a minimum weight fin that would have the planform size needed to provide the needed 
aerodynamic lift and moment authority to control the missile and yet have sufficient stiffness to sat- 
isfy a static divergence speed constraint of 1.15 times the limit speed. Of course, the structure 
would have to be strong enough to withstand the applied loads. 

The second problem addresses the design of the wing for a supersonic transport airplane. 
Candidate planforms would be a clipped delta wing, or a clipped delta wing with cranks in either 
the leading edge or in the trailing edge or in both. Piston theory would be used for this application. 
The wing would be sized to meet some performance requirements, say to develop sufficient lift to 



carry a prescribed payload. A flutter constraint could be applied as well as static deformation con- 
straints. Aerodynamic sensitivity derivatives related to such geometric characteristics as airfoil 
section and planform area could be used. A minimum weight structure would be designed to sat- 
isfy the performance requirements within specified geometric, stress, and flutter constraints. 

CONCLUDING REMARKS 

Two existing aerodynamic theories that readily lend themselves to calculating certain aero- 
dynamic sensitivity derivatives for use in multidisciplinary optimization studies have been reviewed 
briefly. The basic equations relating surface pressure (or lift) to normal wash was given and dis- 
cussed in each case. Although no all-inclusive, rigorous developments of aerodynamic sensitivity 
derivatives were included, the exposition was taken to sufficient depth to provide the foundation 
needed by those wishing to use either of these theories to development aerodynamic sensitivity 
derivatives for use in multidisciplinary optimization studies. In addition, two sample problems are 
suggested in very general terms for studying the process of using aerodynamic sensitivity 
derivatives in optimization studies. 


REFERENCES 

1 . Sobieszczanski-Sobieski, Jaroslaw: The Case for Aerodynamic Sensitivity Analysis. Pro- 
ceeding of Symposium on Sensitivity Analysis in Engineering. NASA CP 2457, pp. 77-96, 
Sep. 25-26, 1986. 

2. Yates, E. Carson, Jr.: Aerodynamic Sensitivities from Subsonic, Sonic, and Supersonic Un- 
steady, Nonplanar Lifting-Surface Theory. NASA TM 100502, Sep. 1987. 

3. Jones Robert T.: Properties of Low- Aspect-Ratio Pointed Wings at Speeds Below and Above 
the Speed of Sound. NACA Rep. 835, 1946. (Supersedes NACA TN 1032.) 

4. Lighthill, M. J.: Oscillating Airfoils at High Mach Number. Journal of the Aeronautical 
Sciences, Vol. 20, No. 6, pp. 402-406, June, 1953. 

5. Ashley, Holt; and Zartarian, Garabed: Piston Theory-A New Aerodynamic Tool for the 
Aeroelastician. Journal of the Aeronautical Sciences, Vol. 23, No. 12, pp. 1 109-1118, Dec. 
1956. 

6. Doggett, Robert. V., Jr.; and Rainey, A. Gerald: Effects of Flexibility on lift and Pitching- 
Moment Characteristics of a Series of low-Aspect-Ratio Wing-Body Combinations at Tran- 
sonic Speeds. NASA TN D-4655, Oct. 1968. (Supersedes declassified NASA TM X-343.) 

7. Rainey, A. Gerald; Hanson, Perry W.; Martin, Dennis J.: Investigation of the Divergence 
Characteristics of Delta-Planform Canard Controls. NASA TR R-235, April 1966. 


15 



8. Morgan, Homer G.; Huckel, Vera; and Runyan, Harry L.: Procedure for Calculating Flutter 
at High Supersonic Speed Including Camber Deflections, and Comparison with Experimental 
Results. NACA TN 4335, Sep. 1958. 

9. Gibbons, M. D.; Soistmann, D. L.; and Bennett, R. M.: Flutter Analysis of Highly Swept 
Delta Wings by Conventional Methods. NASA TM 101530, Nov. 1988. 

10. Heeg, J.; Gilbert, M. G.; and Pototzky, A. S.: Active Control Aerothermoelastic Effects for a 
Conceptual Hypersonic Aircraft. Presented at AIAA Guidance, Navigation and Controls Con- 
ference, AIAA Paper 90-3337, Aug. 1990. 

11. Garrick, I. E.: Some Research on High-Speed Flutter. Third Anglo-American Aeronautical 
Conference (Brighton, England), Royal Aeronautical Society, 1952. 

12. Van Dyke, Milton D.: A Study of Second-Order Supersonic Flow Theory. NACA Rep. 
1081, 1952. (Supersedes NACA TN 2200.) 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0189 


1. AGENCY USE ONLY bUnk) | 2. REPORT OATE j 3. REPORT TYPE AND DATES COVERED 

| August 1991 I Technical Memorandum 

4. TITLE AND SUBTITLE *• NUMBERS 

On Simple Aerodynamic Sensitivity Derivatives for Use in 
Interdisciplinary Optimization 

505-63-50-06 

6. AUTHOR(S) 

Robert V. Doggett, Jr. 


7. PERFORMING ORGANIZATION NAME(S) AND AOORESS(ES) 

NASA Langley Research Center 
Hampton, VA 23665-5225 


PERFORMING ORGANIZATION 
REPORT NUMBER 


9. SPONSORING /MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


10. SPONSORING /MONITORING 
AGENCY REPORT NUMBER 


NASA TM-1041 45 



12a. DISTRIBUTION /AVAILABILITY STATEMENT 

Unclassified-Unlimited 


12b. DISTRIBUTION CODE 


Subject Category 02 


13. ABSTRACT (Maximum 200 words) 

Low-aspect-ratio and piston aerodynamic theories are reviewed as to their use in 
developing aerodynamic sensitivity derivatives for use in multidisciplinary 
optimization applications. The basic equations relating surface pressure (or lift 
and moment) to normal wash are given and discussed briefly for each theory. The 
general means for determining selected sensitivity derivatives are pointed out. In 
addition, some suggestions in very general terms are included as to sample problems 
for use in studying the process of using aerodynamic sensitivity derivatives in 
optimization studies. 


IS. NUMBER OF PAGES 

17 


IB. PRICE CODE 

A0 3 


H i SECURITY CLASSIFICATION 1 18 SECURITY CLASSIFICATION I 19. SECURITY CLASSIFICATION I 20. LIMITATION OF ABSTRACT 

OF RE FORT | OF THIS PAGE OF ABSTRACT I 


14 SUBJECT TERMS 

Aerodynamic sensitivities; Low-aspect-ratio theory; Piston 
theory; Optimization 


j Unclassified 

NSN 7540-01 -280-5S00 


Unclassified 


Unclassified 


PRECEDING PAGE BLANK NOT FILMED 


Standard Form 298 (Rav. 2-89) 

Presented by ANSI Sid 
2M-102 








